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ABSTRACT 



We present the first detection of carbon monoxide (CO) in a damped Lyman-a system (DLA) at z a b s =2.41837 
toward SDSS J143912. 04+111740. 5. We also detected H2 and HD molecules. The measured total column densities (in 
log units) of H 1, H 2 , and CO are 20.10i0.10, 19.38±0.10, and 13.89±0.02, respectively. The molecular fraction, / = 
2A(H 2 )/(iV(H l)+2A(H 2 )) = 0.27+°;^ , is the highest among all known DLAs. The abundances relative to solar of S, Zn, 
Si, and Fe are —0.03+0.12, +0. 16T0.11, —0.86+0.11, and — 1.32±0.11, respectively, indicating a high metal enrichment 
and a depletion pattern onto dust-grains similar to the cold ISM of our Galaxy. The measured JV(CO)/JV(H2) = 3x 10~ 6 
is much less than the conventional CO /H 2 ratio used to convert the CO emission into gaseous mass but is consistent with 
what is measured along translucent sightlines in the Galaxy. The CO rotational excitation temperatures are higher than 
those measured in our Galactic ISM for similar kinetic temperature and density. Using the C I fine structure absorption 
lines, we show that this is a consequence of the excitation being dominated by radiative pumping by the cosmic 
microwave background radiation (CMBR). From the CO excitation temperatures, we derive Tqmbr = 9.15±0.72 K, 
while 9.315±0.007 K is expected from the hot big-bang theory. This is the most precise high-redshift measurement of 
Tcmbr and the first confirmation of the theory using molecular transitions at high redshift. 

Key words. Galaxies: abundances - Quasars: absorption lines - Quasars: individual: SDSS J143912. 04+111740. 5 



1. Introduction 

Damped Lyman-a systems (DLAs) in QSO spectra are 
characterized by very high H 1 column densities, iV(H l)<; 
10 20 cm~ 2 . The inferred metallicities relative to solar vary 
between [Zn/H]= —2.0 and for 2 < z a b s < 3 (e.g. Pettini 
et al. 1997; Prochaska & Wolfe 2002). Therefore DLAs are 
believed to be located in the close vicinity of star-forming 
regions. The dust content in a typical DLA is less than or 
equal to 10% of what is seen in the Galactic ISM for simi- 
lar iV(H 1), however sufficient for favoring the formation of 
H 2 (Ledoux et al. 2003). 

The abundance of H2, the relative populations of the 
H2 rotational levels, and the fine-structure levels of the C 1 
ground-state are used to derive the physical conditions in 
the gas, such as temperature, gas pressure, and ambient ra- 
diation field (Savage et al. 1977; Black & van Dishoeck 1987; 
Jenkins & Tripp 2000; Tumlinson et al. 2002). These con- 
ditions are believed to be driven by the injection of energy 
and momentum through various dynamical and radiative 
processes associated with star formation activity. Thus, de- 
tecting H2 in DLAs at high redshifts is an important step 
forward in understanding the evolution of normal galaxies. 
Detecting other molecules would pioneer interstellar chem- 



istry studies at high redshift (see e.g. Wiklind & Combes 
1995 for the dense ISM component). 



Send offprint requests to: R. Srianand, anand@iucaa.ernet.in 

* Based on observations carried out at the European Southern 
Observatory (ESO), under programme 278.A-5062 with the 
UVES echelle spectrograph installed at the ESO Very Large 
Telescope (VLT), unit Kueyen, on Mount Paranal in Chile. 



In the course of our recently completed Very Large 
Telescope survey for H2 in DLAs, we gathered a sample 
of 13 H2 absorption systems at 1.8 < z < 4.3 out of a 
total of 77 DLAs (Ledoux et al. 2003 & 2006; Petitjean 
et al. 2006; Noterdaeme et al. 2008). Absorption lines of 
HD are detected in one of the DLAs (Varshalovich et al. 
2001) and none show detectable CO absorption. We no- 
ticed a strong preference for HVbearing DLAs being asso- 
ciated with C 1 absorption (Srianand et al. 2005) and having 
high metallicities and large depletion factors (Petitjean et 
al. 2006; Noterdaeme et al. 2008). In the Sloan Digital Sky 
Survey data base, we identified a most promising candi- 
date at z ahs = 2.4185 towards SDSS J143912.04+111740.5 
showing such characteristics. We were allocated 8 hours of 
director discretionary time on the Ultraviolet and Visual 
Echelle Spectrograph (UVES) at the VLT of the European 
Southern Observatory (ESO) to search for CO in addition 
to H2. This observation resulted in the detection of CO 
UV absorption lines that have been elusive for more than a 
quarter century (Varshalovich & Levshakov 1981; Srianand 
& Petitjean 1998; Cui et al. 2005). We also detect H 2 and 
HD absorption lines. In this letter we focus our attention 
on the CO excitation. A detailed analysis of the HD ab- 
sorptions will be presented elsewhere. 
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Fig. 1. A sample of molecular and heavy element absorp- 
tion lines associated with the damped Lyman-a system to- 
ward SDSS J143912. 04+111740.5. The normalized flux is 
given on a velocity scale with the origin at z a b s = 2.41837. 
Smooth curves in each panel show our Voigt profile fits to 
the data. The dashed profiles shown for J< 1 H2 lines are 
the la ranges. Tick marks in each panel indicate the loca- 
tions of the 6 components that are used to fit the J = 3 H2 
absorption lines. The symbol "CO" marks the locations of 
two components detected in S 1 and CO. 



2. Observations 

Both UVES spectrographic arms were used with standard 
dichroic settings and central wavelengths of 390 nm and 
580 nm (or 610 nm) for the blue and red arms, respec- 
tively. The resulting wavelength coverage was 330—710 nm 
with a small gap between 452 and 478 nm. The CCD pixels 
were binned 2x2 and the slit width adjusted to 1", yield- 
ing a resolving power of R = 45000 under seeing conditions 
of ~ 0'.'9. The total exposure time on source exceeded 8 h. 
The data were reduced using the UVES pipeline version 
3.3.1 based on the ESO common pipeline library system 
(M0ller Larsen et al. 2007). Wavelengths were rebinned to 
the heliocentric rest frame and individual scientific expo- 
sures co-added. In the following, we adopt the solar refer- 
ence abundances from Morton (2003). 



3. Analysis 

An asymmetric Lyman-a absorption line, together with 
the corresponding Lyman-/? and 7 absorption lines, indi- 
cates the presence of multiple components. The simultane- 
ous fit to the three H 1 absorption features gives log N(H 1) 
(cm" 2 ) = 17.8, 19.25, 19.20, 19.40, and 20.10±0.10 at the 
velocities of v = -892, -594, -382, -117, and km s _1 
relative to z a b s = 2.41837. The component with maximum 
H 1 column density (at v = Q km s _1 ) also exhibits H2, 
HD, and CO absorption lines spread over ~50 km s _1 (see 
Fig. [TJ. Transitions from the J = 3 H2 rotational level are 
detected in six distinct components. Both HD and CO are 
detected in three and two of these components, respectively. 
Transitions from J — and 1 H2 rotational levels are highly 
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Fig. 2. Voigt profile fits to 12 CO A-X bands detected at 
z abs = 2.4185 towards SDSS J143912.04+111740.5. The 
vertical lines mark the locations of CO absorptions from 
different J levels at the redshift of the main component. 
The profiles in dashed and dotted lines are, respectively, 
for the first and second components (at z a b s = 2.41837 and 
2.41847) and the total profile is shown in solid line. 



saturated, but accurate integrated column densities can be 
derived from damped wings. 

Absorption lines from N I, O I, C I, C I*, C I**, Mg I, 
Ar I, S I, S II, Si II, Fe n, Zn II, Al II, and Ni n are seen spread 
over up to 950 km s _1 . The total gaseous abundances rela- 
tive to solar are -0.03±0.12, +0.16±0.11, -0.86±0.11, and 
— 1.32±0.11 for S, Zn, Si, and Fe, respectively The abun- 
dances of S and Zn are consistent with the gas abundance 
being close to solar, while the relative depletions of Si and 
Fe are similar to those in cold gas in the diffuse Galactic 
ISM. Absorption lines from the three C I ground-state fine- 
structure levels are detected in numerous transitions in the 
main H I component. 

In Fig. [1] we show a few transitions from H2 (J = 0, 
1 and 3), S 1, C 1, HD and CO. By fitting the damping 
wings of J = and 1 transitions, we measured log V(H2, 
J=0) = 18.90±0.10 and log N(R 2 , J=l) = 19.18±0.10. We 
derived an excitation temperature of T01 = 105^32 K, which 
is usually a good indicator of the average kinetic temper- 
ature of the gas. The mean molecular fraction of the gas 
is / = 2V(H 2 )/[2V(H 2 )+V(H 1)] = 0.27±g;pg. This is the 
highest value measured to date in a high- 2: DLA. 



3.1. The CO molecules and the CO/H2 ratio at high z 

Carbon monoxide absorption was detected in several bands 
(see Fig. [5]). We used the four bands that are redshifted 
outside the Lyman-a forest to derive CO column densities. 
Absorptions from different CO rotational levels are located 
close to each other, but the resolution of the data is high 
enough to dcblcnd them. A single-component Voigt profile 
fit gives z a bs = 2.41837 for the main CO component. This is 
coincident within error (<1 km s _1 ) with the strongest S I 
component (See Fig. [2]). A weaker S 1 satellite component is 
present at z a b s = 2.41853. We therefore added a second CO 
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Fig. 3. The CO excitation diagram. A straight line with 
slope 1/(T GX In 10) indicates thermalization of the lev- 
els. The diagram is given for the main CO component at 
Znbs — 2.41837. The three lines give the mean and la range 
obtained from Tqi, T02, and T±2. The diagram is compatible 
with thermalization by a black-body radiation of temper- 
ature 9.15±0.72 K when T C mbr = 9.315±0.007 K (long 
dashed line) is expected at z a ^ s = 2.4185 from the hot big- 
bang theory. 



component with a redshift fixed to the value of the second 
S 1 component. Doppler parameters were left free to vary 
and the best fit was obtained for t>=1.5 km s" 1 , but the col- 
umn density value in the main component is not very sensi- 
tive to the exact value for b > 0.5 km s" 1 . Results are shown 
in Fig. H Column densities are log iV(CO) = 13.27±0.03, 
13.48±0.02, and 13.18±0.06, respectively, for J = 0, 1 and 
2 in the main component. Column densities for J — and 1 
are 12.75±0.05 and 12.89±0.08 respectively for the second 
component when b = 1.5 km s" 1 . 

We measured N{CO)/N(R 2 ) = 3xl0" 6 . This is sim- 
ilar to or slightly higher than what is measured along 
the Galactic sightlines with similar molecular fraction and 
along Galactic sightlines with similar N(H.2 ) (see Figs. 4 
and 5 of Burgh et al. 2007). This is much less than the 
CO/H2 ratio of about 10" 4 derived for dense molecular 
clouds (e.g Lacy et al. 1994). This strongly suggests that 
the physical conditions in the gas are similar to those in 
the diffuse Galactic ISM. 

3,2. CO rotational excitation 

The excitation temperatures derived from the population 
ratios of the different rotational levels are Tqi = 9.11 ± 1.23 
K, T12 = 9.19±1.21K, andT 02 = 9.16±0.77 K for the main 
component where the errors come from the fitting uncer- 
tainties. Additional rms deviations come from uncertainties 
in the continuum placement and from the allowed range for 
the Doppler parameter of the second component. We esti- 
mate these to be ^0.21, ^0.37, and ^0.18 K around the 
three excitation temperatures. The populations of the three 
rotational levels are thus consistent with a single excitation 
temperature, T ox = 9.15 ± 0.72 K (see Fig. 3), suggesting 
that a single mechanism controls the level populations. 
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Fig. 4. Comparison of rotational excitation temperatures 
of CO and H2 . Filled circles are for the measurements 
from Galactic diffuse ISM (Burgh et al. 2007). The filled 
triangle is for our measurement at z — 2.4183 towards 
SDSS J143912.04+111740.5. 



This excitation temperature is more than twice higher 
than the excitation temperature measured in the diffuse 
Galactic ISM (see Fig. @|, at similar H 2 T m (Burgh et al. 
2007). Since we have seen that the physical conditions in 
the DLA are similar than in the diffuse Galactic ISM, we 
can suspect that the process responsible for this high ex- 
citation temperature is specific to the high redshift of the 
system. Interestingly, the population ratios in the second 
component are also consistent with this value of T cx , albeit 
with large errors. 

The presence of absorptions from the three fine- 
structure levels of the C 1 ground state allows us to probe 
the physical state of the gas further (Srianand et al. 2000). 
By simultaneous Voigt profile fitting of the absorption 
lines, we derived log N(C 1) = 14.26±0.01, log N(C I*) = 
14.02±0.02, and log N{C I**) = 13.10±0.02 for the velocity 
component that corresponds to the strongest CO compo- 
nent. Assuming no contribution from the CMBR, we obtain 
a conservative hydrogen density range of 87-135 cm" 3 and 
52-84 cm" 3 from the population ratios N(C i*)/N(C 1) 
and N(C l**)/JV(G 1), respectively. A more realistic range, 
45-62 cm" 3 , is obtained if we use appropriate excitation 
by the CMBR with a temperature expected from the hot 
big-bang theory. From the observed / value, we derived a 
conservative (< 25 cm -3 ) and a realistic (< 12 cm" 3 ) value 
for the H 2 density (tih 2 = /ijhi/ (2 — /))• The above values 
are strict upper limits, because we have ignored the contri- 
bution of UV photons from stars in this galaxy to the C 1 
excitation. 

We ran the statistical equilibrium radiative transfer 
code RADEX, available on line (van der Tak et al. 2007), 
and found that for the kinetic temperature, T = 105 K, de- 
rived from H2, the collisional contributions to Tqi and T12 
are <5% and < 2% for uh 2 < 25 cm -3 . The corresponding- 
values are < 3% and <1% for nn 2 < 12 cm" 3 . Thus the 
collisional excitation of CO by H2 is negligible. Collisions 
with H contribute little to the CO excitation compared 
to collisions with H2 in astrophysical conditions (Green & 
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Fig. 5. Measurements of T(CMBR) at various z. The star 
with errorbars is for the measurement based on CO pre- 
sented here. Our earlier measurement using fine-structure 
lines of neutral carbon, 6.0 < Tcmbr < 14.0 K, at z = 
2.33771 is indicated by the long vertical bar (Srianand et 
al. 2000). Measurement by Molaro et al. (2002) is marked by 
an open circle. Upper limits are measurements using C I or 
C II* from the literature (squares) and using C I from our 
UVES sample (hexagons). The dashed line is the predic- 
tion from the hot big-bang theory, Tcmbk(z)=Tcubr(z = 
0) x (1 + z). Tcmbr measurement at z = is based on the 
COBE determination (Mather et al. 1999). 



Thaddeus, 1976; Shepler et al. 2007). This means that the 
CO excitation is dominated by CMBR, so we conclude that 
T ex = Tcmbr = 9.15±0.72K. 



4. Conclusion 

The CMBR is an important source of excitation for those 
species with transitions in the sub-millimeter range. This is 
the case for atomic species whose ground state splits into 
several fine-structure levels and of molecules that can be 
excited in their rotational levels. If the relative level popu- 
lations are thermalized by the CMBR, then the excitation 
temperature gives the temperature of the black-body ra- 
diation. It has long been proposed to measure the relative 
populations of such atomic levels in quasar absorption lines 
to derive T CM br at high redshift (Bahcall & Wolf 1968). 
In Fig. [5]we combine our precise measurement of Tcmbr, 
the 51 new upper limits obtained using C I and C I* ab- 
sorption lines detected towards QSOs in our UVES sam- 
ple (Srianand et al. 2005; Noterdaeme et al. 2007a, 2007b; 
Ledoux et al. 2006), and measurements reported from the 
literature (Meyer et al. 1986; Songaila et al. 1994; Lu et al. 
1996; Ge et al. 1997; Roth & Bauer, 1999; Molaro et al. 
2002, Cui et al. 2005). Upper limits are obtained assuming 
CMBR as the only source of excitation. Our precise mea- 
surements using CO and the new upper limits using C I are 
consistent with the adiabatic evolution of Tcmbr expected 
in the standard big-bang model (Fig. [5]). 

The CN molecule has proven to be a remarkable ther- 
mometer of the CMBR in our Galaxy. It has been used 
for precise measurement of Tcmbr. in different directions 



(Meyer et al. 1985; Kaiser et al. 1990). Wiklind & Combes 
(1996) obtained Tcmbr < 6 K at z = 0.885 using the ab- 
sorption lines of CS, H 13 CO + , and N2H+. Carbon monox- 
ide in diffuse gas provides an interesting possibility for 
measurements at high redshift, as the rotational energies 
between different rotational levels are close to Tcmbr at 
z > 1. 

Following careful selection of the target, based on in- 
tensive observations at ESO-VLT, we have made the first 
detection of carbon monoxyde molecules in the diffuse ISM 
at high redshift. The analysis presented here pioneers inter- 
stellar chemistry studies at high redshift and demonstrates 
that, together with the detection of other molecules such 
as HD or CH, it will be possible to tackle important cos- 
mological issues. 
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